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STATIONARY MAX-STABLE FIELDS ASSOCIATED TO 
NEGATIVE DEFINITE FUNCTIONS 

By Zakhar Kabluchko, Martin Schlather and Laurens de Haan 

University of Gottingen, University of Gottingen and 
Erasmus University Rotterdam 

Let Wi,i € N, be independent copies of a zero- mean Gaussian 
process {W(t),t £ K d } with stationary increments and variance o~ 2 (t). 
Independently of Wi, let ^-=i De a P°i ss °n point process on the 
real line with intensity e~ y dy. We show that the law of the random 
family of functions {Vi(-),i G N}, where Vi(t) = Ui + Wi(t) - a 2 (t)/2, 
is translation invariant. In particular, the process rj(t) = \/ € *l 1 Vi(t) 
is a stationary max-stable process with standard Gumbel margins. 
The process r\ arises as a limit of a suitably normalized and rescaled 
pointwise maximum of n i.i.d. stationary Gaussian processes as n — > 
oo if and only if W is a (nonisotropic) fractional Brownian motion 
on R d . Under suitable conditions on W, the process rj has a mixed 
moving maxima representation. 

1. Introduction. A stochastic process {rj(t),t £ R d } is called max-stable 
if, for any n £ N, the process {Vfc=i t £ nas the same law as {r](t) + 
logn,t G W 1 }, where r)i,...,r] n are independent copies of r\. It follows from 
this definition that the marginal distributions of i] are of the form P [??(£) < 
x] = exp(— e - x + b i t )^j and, more generally, the finite-dimensional distributions 
of r\ are multivariate max-stable distributions of Gumbel type [26]. Max- 
stable processes have been studied in [8, 10, 12, 16, 29] and [9], Part III. Note 
that it is common to consider max-stable processes with Frechet (rather than 
Gumbel) marginals, so most authors work with the process e 11 instead of i]. 

A general description of stationary max-stable processes in terms of non- 
singular flows on measure spaces was given in [12]. A usual approach to 
constructing examples of such processes is to use some sort of moving max- 
ima (or, more generally, mixed moving maxima) representation; sec [11, 14, 
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27, 31]. Another family of examples, based on stationary random processes, 
was given in [27]. Contrary to the mixed moving maxima processes, which 
were shown to be mixing in [30], this family entails a nonvanishing large- 
distance dependence within the max-stable process. 

In this paper, we are mainly interested in a remarkable stationary max- 
stable process constructed by Brown and Resnick in [4] . Let us recall part of 
their result (see also Section 9.8 in [9] for the two-sided version given here). 

Theorem 1. Let W%,i G N, be independent copies of a standard Brown- 
ian motion {W(t),t G K} and, independently o/Wj, letJ2iZi °~Ui be a Poisson 
point process on R with intensity e~ v dy. Then, the process 

oo 

(1) r,(t) = \/(U i + W i (t)-\t\/2) 

i=l 

is a stationary max-stable process with standard Gumbel margins. 

A natural question arises as to whether further stationary max-stable 
processes can be constructed by replacing, in the above construction, the 
drifted Brownian motion Wit) — \t\/2 by other stochastic processes. Thus, 
we are interested in stochastic processes {£(£), i G M. d } having the property 
that the process rj(t) = V~i(E^i + is stationary, where Ui,i G N, are 

as above and G N, are independent copies of £. We call such processes 
£ Brown-Resnick stationary] see Section 2 for a more precise definition. In 
[4], two different proofs of Theorem 1 were given. One of them is based 
on the fact that e~ y dy is an invariant measure for the Brownian motion 
with drift — 1/2 and can be extended to show that some classes of processes 
with Markov property are Brown-Resnick stationary; see [5, 30]. The other 
proof, which uses the connection with the extreme- value theory of Gaussian 
processes, will be discussed later in Sections 6 and 8. 

We are going to show that any Gaussian process with stationary incre- 
ments becomes Brown-Resnick stationary after subtracting an appropriate 
drift term. Recall that a random process {W(t),t G ffi d } is said to have 
stationary increments if the law of {W(t + to) — W(to), t G U. d } does not de- 
pend on the choice of to G M d . If W is a Gaussian process with stationary 
increments (always supposed to have zero mean), then its law is completely 
characterized by what we shall call the variogram 

-f(t)=EiWit + t )-W(t )) 2 , t£R d , 

and the variance a 2 it) = Var Wit). It is well known that a function 7 : M d — > 
[0, 00) with 7(0) = is a variogram of some Gaussian process with stationary 
increments if and only if it is negative definite. The latter condition means 
that "f{—t) = 7(t) for every t G M d and Z)ij=i a i a j"f{U ~ *i) — for every 
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ti, . . . , t n E M. d and a\, . . . , a n E R satisfying X)it=i a « = 0> see [2] for more on 
negative definite functions. Examples of Gaussian processes with stationary 
increments are provided by, for example, stationary Gaussian processes, their 
integrals (if d = 1) and fractional (Levy) Brownian motions, the latter being 
characterized by W(0) = and j(t) = \\t\\ a for some a £ (0,2]. Here, ||t|| 
denotes the Euclidean norm of the vector t. 

Theorem 2. Let Wi,i E N, be independent copies of a Gaussian process 
{W(t),t E M. d } with stationary increments, variance o~ 2 {t) and variogram 
j(t). Independently ofWi, let J2i^i &Ui be a Poisson point process on R with 
intensity e~ y dy . The process 

oo 

(2) n (t ) = \J(jj i + Wi{ t)-a 2 {t)/2) 

i=\ 

is then a stationary max-stable process with standard Gumbel margins. The 
law of n depends only on the variogram -y. 

The process r\ defined above will be called the Brown-Resnick process as- 
sociated to the variogram 7. If j(t) = \t\ [i.e., if the corresponding Gaussian 
process W, under W(0) = 0, is a standard Brownian motion], then we recover 
the process of Theorem 1, originally considered in [4]. The Brown-Resnick 
process corresponding to the variogram "f(ti, . . . ,tj) = \t\\ + • • • + |i<j| was 
used as a model of extreme spatial rainfall in [6] and [13]. Another natural 
class of random processes, having the advantage of being isotropic, can be 
obtained by taking j(t) = \\t\\ a , t E M. d , for some a E (0,2]. If a = 2, the cor- 
responding drifted Gaussian process W(t) — a 2 (t)/2 is a "random parabola" 
of the form W(t) = (t,N) — ||t|| 2 /2, where the random vector iV has the 
standard Gaussian distribution on Mr and we recover a process introduced 
in [15] and [17]; see also [11]. If 7 is bounded, then the process W can be 
chosen to be stationary (after changing the variance and without changing 
the variogram; see, e.g., Proposition 7.13 in [2]) and r\ belongs to the class 
of max-stable processes considered in Theorem 2 of [27]. 

Different Gaussian processes with stationary increments may have the 
same variogram. For example, let {W(t),t E M} be a standard Brownian 
motion and let / E L 2 (R). The process W f (t)=W(t) + J R f(s)dW(s) then 
has the same variogram j(t) = \t\ as W and it is not difficult to see that the 
laws of Wf and W g coincide if and only if / = g a.s. The fact that different 
processes with the same variogram lead to the same 77 is quite surprising, 
even in the particular case mentioned above. 

The Brown-Resnick processes defined in Theorem 2 have no a priori con- 
nection to mixed moving maxima processes mentioned at the beginning of 
the paper. It was asked in [30] if the original Brown-Resnick process cor- 
responding to j(t) = \t\ has a representation as a mixed moving maxima 
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process. We shall show in Section 5 that the answer is affirmative. More 
generally, it will be shown that the Brown-Resnick process corresponding 
to a Gaussian process W with stationary increments has a mixed moving 
maxima representation provided that lim|| t ||_ >00 (W / (t) — a 2 (t)/2) = —oo a.s. 

The paper is organized as follows. In Section 2, we introduce the notion 
of Brown-Resnick stationarity. In Section 3, we prove a general criterion 
which allows one to decide whether a given random process £ has the prop- 
erty of Brown-Resnick stationarity in terms of the Laplace transform of the 
finite-dimensional distributions of £. This criterion is then used in Section 
4 to prove Theorem 2. In Section 5, we show that Brown-Resnick processes 
of Theorem 2 have a mixed moving maxima representation under some con- 
ditions on the variogram 7. In Sections 6 and 7, we study, generalizing [4], 
extremes of a large number of independent Gaussian processes. An alterna- 
tive proof of Theorem 2, in the case W(0) = 0, is given in Section 8. 

Remark 3. Two objects will appear frequently in our considerations: 
the Poisson point process Ya=i &Ui with intensity e~ y dy on M and the stan- 
dard Gumbel distribution exp(—e~ y ), which is the distribution of max^g^ Ui. 
The transformation y 1— ► e y allows us to switch from Gumbel to the more 
common Frechet notation. That is, if Y is a random variable with standard 
Gumbel distribution, then Z = e Y has standard Frechet distribution, mean- 
ing that ¥[Z < z]= exp(— 1/z), z > 0. Further, J2i=i $ e u i ^ s a Poisson point 
process on (0, 00) with intensity dz/z 2 . Thus, if 77 is a max-stable process, 
as defined at the beginning of the paper, then the process e 11 is max-stable 
in the usual sense [8]. 

2. Brown Resnick stationarity property. Let € N, be independent 
copies of a random process {£(t),t G M. d } satisfying 

(3) Ee m < 00 for all t G R d . 

Further, let $Ui be a Poisson point process on M with intensity e~ y dy, 
independent of the family i € N. Define a process {rj(t),t G M. d } by 

00 

(4) ri(t) = y(Ui + ti(t))- 

i=i 

The process 77 is necessarily max-stable [8]. To give a short proof of this 
fact, let 771, . . . ,77„ be independent copies of 77, constructed by starting with 
k i k = 1, . . . ,n, and i G N, k = 1, . . . ,n, all objects being inde- 
pendent. The superposition J2k=i J2t=i ^u, k is then a Poisson point process 
on R with intensity ne~ y dy = e -(y~^ n ) dy. Hence, £fc=i S Ui k _i og „ has 
the law of the Poisson point process with intensity e~ y dy. So, the process 
Vfc=i % — logn has the same law as 77, which proves the max-stability of 77. 
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By [8], the converse is also true: any stochastically continuous max-stable 
process r\ is of the form (4) for some process £. The finite-dimensional dis- 
tributions of r] were computed in [8]: given ti, . . . , t n G R d and y\, . . . , y n G R, 
we have 

(5) PfoOi) <yi,---,r)(t n ) <y n ] = exp{-Eexp max (£(U) - yi)\ . 

y. 1=1, ...,n ) 

In particular, condition (3) ensures that for every t G R d , rj(t) is finite a.s. 
We are interested in processes £ leading to a stationary process rj. 

Definition 4. A stochastic process {£(t),t G R d } satisfying (3) is called 
Brown-Resnick stationary if the process rj defined in (4) is stationary. 

It is trivial that every stationary process satisfying (3) is Brown-Resnick 
stationary. However, the converse is not true: by a result of [4], the non- 
stationary process £(i) = W(t) — \t\/2, where {W(t),t G R} is a standard 
Brownian motion, is Brown-Resnick stationary. The next proposition gives 
an equivalent, but perhaps more natural, version of Definition 4. 

Proposition 5. A process {£(t),t G R d } which satisfies (3) is Brown- 
Resnick stationary if and only if ^t/i ^ s a translation invariant 
Poisson point process on the space E = M Rd . 

Before we can start the proof, we need to introduce some notation. We 
endow E = R* , the space of real- valued functions on M. d , with the product 
ex-algebra 13(E) generated by the finite-dimensional cylinder sets, that is, by 
the sets of the form 

(6) C tlt ..., tn (B) = {/ : R d -> R : (f(h), f(t n )) G B}, 

where t\ , . . . , t n G R d and B is a Borel set in R n . If the processes £j have 
continuous sample paths, then E = C(R d ), the space of continuous func- 
tions, could be considered as well. Let A4(E) be the space of all measures 
on E which have the form fi = $fi f° r some G E and which are lo- 
cally finite [i.e., finite on all cylinder sets of the form (6) with bounded 
B]. We endow M.(E) with the a-algebra B(M(E)) generated by the maps 
F tl ,...,t n ;B-M(E) ^N U{c»}, /jh ii(C tu ,„ u (B)). A point process on E 
is a random variable O : O — > M (E) , defined on some probability space f2 
and taking values in M(E). Also, recall (see [19, 26]) that for a locally fi- 
nite measure A on E, a Poisson point process with intensity A is a point 
process Q-.n^M(E) such that @(-)(A) ~ Poiss(A(A)) for each AeB(E), 
A(A) < 00, and the random variables @(-)(Ai), i G N, are independent pro- 
vided Ai G B(E) are disjoint. 
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We define a family of operators Th'.M.(E) —* M(E), h G R d , as follows: 
for fi = YliliSfi &M(E), we define T h (fi) = Y^Li $fi(-+h)- A P oint process 
on E is called translation invariant if its distribution, viewed as a probability 
measure on M(E), is invariant with respect to the family T^. A measure 
A on the space E is called translation invariant if, for every A G B{E) and 
every h G R d , we have A(A) = A({/(- + h) :/ G A}). A Poisson point process 
on i? is translation invariant if and only if its intensity measure A is 
translation invariant. 

Proof of Proposition 5. Let be the law of £ on the space E = R Rd . 
Define a map it : R x E -> £ by 7r(J7, £(•)) = [/ + £(•) and let A be the push- 
forward of the measure e~ y dy x dRt by the map 7r [i.e., for ^4 G B(E), define 
A.(A) = j 7r ^ 1 ^e~ y dy x dRt]. We show that condition (3) implies that the 

measure A is locally finite. To this end, take t G R d and let A^k = {/ G 
E:f{t)>k}, keZ. Then, 

A(A tjk ) = f e~"P[e(t) >k-y]dy = e~ k [ e*P[f (t) > z] dz, 
Jm. Jm 

which is finite, by (3). Since any bounded cylinder set is contained in some 
At t k, the measure A is locally finite. Since Ufcgz^fc = the measure A is 
cr-finite. 

The random measure Ya^i ma y be viewed as a Poisson point 

process on R x E with intensity e~ y dy x d¥^. Therefore, by a general map- 
ping theorem (see [19]), Yi^i^Ui+d(-) is a Poisson point process on E with 
intensity measure A. Given t±, . . . , t n G R d , yi, . . . ,y n G R and denoting = 
R n \ X" =1 (-oo, yj ], we have 

Pfo(tl) < 2/1, • • -Mtn) < Vn] = n& G N : f/j + &(•) G C tl ,.„ )tB (S)] 

(7) 

= exp(-A(C tl) ...,* n (£))). 

Now, suppose that the point process Ya^=i ^+^(0 * s translation invariant. 
It follows that its intensity measure A is translation invariant. Equation (7) 
then implies that the process r/ is stationary. Conversely, if r\ is stationary, 
then, again using (7), we obtain that 

l+h,...,t n +h 

for every set B of the form W 1 \ X !^ =1 (— oo, g/j] and every h G R d . The 
translation invariance of A follows from this, using uniqueness of extension 
of measures and the <7-finiteness of A. □ 
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3. A general stationarity criterion. In this section, we prove a general 
criterion for the Brown-Resnick stationarity of a given process in terms of 
Laplace transforms of its finite-dimensional distributions. Let {£(t),t G R rf } 
be a random process satisfying (3). For t\,...,t n G R d , denote by ^ti,...,t n 
the distribution of the random vector (£(ii), • ■ • , £(t n )). An application of 
Holder's inequality shows that the Laplace transform of the measure ¥ tl ... tn , 
defined by 

<p tl ,...,u(u 1 ,...,u n )= I e u ^ + - +u " x "dP tl _ tn ( Xl ,...,x n ), 

JR n 

is finite provided U{ G [0, 1], Ya=i ui<1. 

Proposition 6. A random process {£(£), i GR rf } satisfying the moment 
condition (3) is Brown-Resnick stationary if and only if 

(8) ^i,...,tn(«l, ...,Un) = (p tl +h,...,t n +h(ui, • • • , «n) 

for every h, t±, . . . , t n G R d and any u\, . . . , u n G [0, 1] satisfying J2?=i 1H = 1. 
We need the following lemma on the uniqueness of the Laplace transform. 

Lemma 7. Let /ii and ^ be two finite measures on R n with Laplace 
transforms ipi(t) = / R „ e^'^ dfj,\(s) and ip2(t) = Jk« e^' s ^ dfiz^s) such that ipi 
and ip2 are finite and equal on some open set D C R n . Then, [i\ = fi2- 

Proof. If ipi and ip2 are finite on D, then they are finite on the com- 
plexification of D, that is, on the set D c = {t G C n :Rei G D}. Since ip± and 
ijj2 are analytic functions coinciding on D, they must coincide on D c . Let 
to G D. Then, s i— ► ip\(to + is) is the characteristic function of the finite mea- 
sure e^ '''* dfii(-). Now, ipi(to + is) = ^2(^0 + and the fact that a finite 
measure is uniquely determined by its characteristic function together imply 
that e< to >-> = e<*°''> d// 2 (-)- Hence, m = fi 2 - □ 

Proof of Proposition 6. We use the notation of the previous section. 
Our goal is to show that the intensity measure A is translation invariant 
if and only if (8) holds. For a set BcK" and x G R, let B + x = B + 
(x,x, .. . ,x). For a cylinder set Ct lt ...,t n {B) [recall (6)], we have 

A(C tl _ tn (B)) 

= [ e x F tl _ tn (B + x)dx 

JR 

= e x l B +x(yi,---,yn)dF tlt ... ttn (y 1 ,...,y n )dx 

JR JR n 
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e vi e *-vi i B+x _ yi (0,2/2 - yi, ■ •■ ,y n - m) d$'ti,...,t n (yi' ■■■,yn)dx 
e yi e z i B + z (o, y2-yi,---,y n -yi) dF tly .. jtn (yi,- ■ ■ , y n ) dz. 

Consider a measure lMt lt ...,t n on ^ n > defined on Borel sets A C M. n by 

fH u ...,tn( A ) = / e Vl l A (0,y 2 - yi, . . . ,y n - yi) d¥ tl tn (yi, . . . ,y n ). 

Note that, by (3), we have (Hi,...,tn(A) < Ee^* 1 ^ < oo and therefore the mea- 
sure fJrtt,...,t„ is finite. The measure (H u ..,,t„ m ay be viewed as a type of 
exponentially weighted projection of the measure t n ° n t° the (n — 1)- 
dimensional hyperplane {(xi)f =1 G M. n :x± = 0}. We have 

(9) A(C tl _ tn (B))= [ e z fMl _ tn (B + z)dz. 

The Laplace transform of Ht\ ...tn is gi yen by 
ipt 1 ,...,t n (ui,...,u n ) 

e tfi e « 2 (tt,-vi)+-+«»( tf »-yi) dP tl + n ( yi , ...,y n ) 



(10) 

where (fa ...t n is the Laplace transform of the measure t n - Note that 
ifax,...,tn does not depend on u\. 

Now, suppose that (8) holds. We then obtain 

(11) ^ii,...,t„( U l> ■■■,Un)= lf) tl +h,...,t n +h(ui, ...,U n ) 

provided that Uj € [0,1], J2i=2 u i — 1j which, by Lemma 7, implies that 
Mti,...,t„ = fH!+h,...,tn+h an d hence, by (9), 

(12) A(C tl+h ,..., tn+h (B)) = A(C tl ,..., tn (B)). 

This proves the translation invariance of A on the semi-ring of the cylinder 
sets. Using the theorem on the uniqueness of the extension of measures and 
the fact that A is cr-finite, we obtain the translation invariance of A on the 
whole cr-algebra 13(E). 

Now, suppose that A is translation invariant. It follows that (12) holds 
and thus, using (9), 

/ e z fi tl ,...,t n (B + z)dz= e z ' fi tl +h,...,t n +h(B + z) dz 
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for every Borel set B C W 1 and every h,t\, . . . ,t n £ R d . Since the measure 
fMi,...,t„ is concentrated on the hyperplane {(xj)" =1 G R n : x\ = 0}, it follows 
that, actually, Hti,...,t n = ^t-i_+h,...,t n +h- By considering the Laplace trans- 
forms, we obtain that (11) holds, from which (8) follows. This completes 
the proof. □ 

As an immediate consequence of the above proposition, we obtain the 
following nontrivial corollaries: 

Corollary 8 . Let {£' (t) , t e R d } and {£"(£), t £ R d } be two independent 
processes, both having the Brown- Re snick stationarity property. The process 
£' + £" is then also Brown-Resnick stationary. 

Corollary 9. Let {£i(t),t e R dl } and {&(t),t £ R d2 } be independent 
Brown-Resnick stationary processes. The process {£(ii, £2), (ti 5 ^2) G R dl+d-2 } 
defined by £(ii,t2) = £i(ii) + £2(^2) is then Brown-Resnick stationary. 

4. Max-stable processes associated to variograms. 

Theorem 10. Let {W(t),t £ R d } be a Gaussian process with stationary 
increments and variance o~ 2 (t). The process £(t) = W(t) — o~ 2 {t)/2 is then 
Brown-Resnick stationary. 

Proof. Recall our standing assumption E,(W(t)) = 0. It follows from 
the definition of the variogram j(t) = E,(W(t) — W(0)) 2 that we have 

Cov(W(t), W(s)) = o- 2 (t)/2 + a 2 {s)/2 - 7 (t - s)/2. 

We are going to apply Proposition 6 to £(i). Note that Ee^ = 1, which 
shows that (3) is satisfied. We need to prove that (8) holds. The distribu- 
tion Pt lr ..,t n of the random vector (£(£i), . . . , £(in)) is a multivariate Gaus- 
sian distribution whose expectation vector (//i)i=i n an d covariance matrix 
{<Jij)i,j=i,...,n are given, respectively, by 

(13) ^ = -a 2 {ti)/2 and a {j = a 2 (U)/2 + a 2 (tj)/2 - 7^ - tj)/2. 
The Laplace transform of Rt lt ...,t n is given by 

/ n 1 n \ 

(14) (p tu ,„ >tn (ui,...,u n ) =exp X^AW+o a ij u i u j)- 

\i=i 1 i,j=i J 

Let ui, . . . , u n £ [0, 1] satisfy Y^i=i u i = 1- substituting ui = 1 — J27=2 u i 
into (14) and using (13), we obtain that 

(15) <p tlt ...,t n (ui,...,u n )=exp(L + lQ), 
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where L = L ilv .. iin (n 2 , . . . ,u n ) and Q = Qt lt ...,t 7l (w 2 , . . . ,u n ) are the linear 
part and the quadratic part, respectively (the constant term is easily seen 
to be zero). The linear part is given by 

n 1 n 

(16) L = ^(Mi - Ml + °ii ~ = ~ n ~ 

i=2 i=2 

The quadratic part is easily seen to be 
n 

Q= (°X? ~ ° U ~ ai 3 + a ll) u i u j 

= 2 E (7(<i - <l) + 7(*j - <l) - 7(*j - U))uiUj. 

Thus, both terms L and Q do not change if one replaces ti,...,t n by t\ + 
h, . . . ,t n + h. Hence, (8) holds and the proof is complete. □ 

Proposition 11. Let W and W" be two Gaussian processes with sta- 
tionary increments, having the same variogram ^{t) and possibly different 
variances cr' 2 (t) and a" 2 (t). Let A' (resp., A") be the intensity of the Pois- 
son point process constructed as in Section 2 with £ replaced by W — a' 2 /2 
(resp., W"-a" 2 /2). Then, A' = A" . 

Proof. Formulas (15), (16) and (17) of the previous proof show that 
</>ti,...,t„ =¥>ti, ...,t„> which, by (10), implies that ^ lv .., in = ^t 1 ,...,t n - Here > a11 
objects marked with ' (resp., ") correspond to W (resp., W"). Lemma 7 
yields fi' tl t = fj.'^ t . Now, (9) shows that for every cylinder set Ct lt ...j n (B), 
we have 

A\C tl ,..., tn (B)) = A"(C tl ,... itn (B)). 
To finish the proof, use the a-finiteness of A' and A". □ 

Remark 12. Given a Gaussian process W with stationary increments, 
it will often be convenient to replace it by the process W(t) = W(t) — W(0) 
having the same variogram 7 as W and W'(O) = 0. Note that the variance 
of the process W is given by a 2 (t) = j(t). 

Proof of Theorem 2. The stationarity of r\ follows from Theorem 10, 
whereas the max-stability was proven in the discussion following (4). The 
fact that r](t) is standard Gumbel for each t € M. d follows from (5). Finally, 
the last claim of the theorem follows from Proposition 11. □ 
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Proposition 13. If all Gaussian processes in Theorem 2 have contin- 
uous sample paths, then the process ij is also sample- continuous. 

PROOF. Let K C R d be bounded. We use the notation £(t) = W(t) - 
a 2 {t)/2 and = Wi(t) - cr 2 (t)/2. First, we show that for every k £ Z, the 
random set 

I k = \ieN:sup(Ui + £i(t))>k} 

is a.s. finite. Indeed, the cardinality of I\. is Poisson distributed with some 
(maybe infinite) intensity We have 



/OO p -, /"OO r 

e~ z ¥ sup^(t) >k-zdz<l + e z F sup^(t) > jfe + z 
-oo HeK J JO l t&K 



dz. 



Since the process £ is Gaussian with continuous paths, a result of [20] (or 
see Corollary 3.2 of [22]) states that Eexp{e(sup t6 ^ £(i)) 2 } < oo for some 
small e > 0. Hence, < oo and, consequently, Ik is finite a.s. 

We now show that r\ is continuous a.s. Let Ak, k £ Z, be the random event 
infi g /<(Lr + £i(i)) > fc. Note that IP[Ufcez^fc] = 1- If; sa Y> ^fe occurs, then 

7/(i)= \/ ^K. 

iG/ fe U{l} 

It follows that 77, being a pointwise maximum of a finite number of contin- 
uous functions, is itself continuous. □ 

5. Representation as mixed moving maxima process. We are now going 
to show that under some conditions on the underlying variogram 7, the 
Brown-Resnick process 77 has a representation as a mixed moving maxima 
process. First, we recall a definition of mixed moving maxima processes as 
given in [27, 30]; see also [14, 29, 31]. Let {F(t),t £ R d } be a measurable 
process and suppose that E J Rd e F ^dt < 00. Let Y^i^(ti,yi) De a Poisson 
point process on M. d x R with intensity e~ y dt dy (dt is the Lebesgue measure 
on R d ) and let F^i E N, be independent copies of F. A process of the form 



r ] {t) = \J(F i {t-t i )+y i ), t£R d , 
i=l 

is called a mixed moving maxima process. It is convenient to think of F{ as a 
random mark attached to the point (ti,yi). The process r\ is stationary and 
max-stable; its finite-dimensional distributions are given by 

P[??( s i) ^ z \i ■ ■ ■ :V( s n) < z n ] = exp< — E / exp max (F(sj — t) — Zj) dt >, 

I m d i=i,-,« J 

where si, . . . , s n € R d , zi, . . . , ^ n € R and E denotes the expectation with 
respect to the law of F (see, e.g., [29]). 
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Theorem 14. Let {W(t),t G R } be a sample- continuous Gaussian pro- 
cess with stationary increments and variance o~ 2 {t). Suppose that 

(18) lim (W(t) -a 2 (t)/2) = -oo a.s. 

1 1 1 1 1 — >oo 

The Brown-Resnick process r\ defined in Theorem 2 then has a representa- 
tion as a mixed moving maxima process. 

Proof. Recall that YA^i^Ui 1S a Poisson point process on R with in- 
tensity e~ y dy and Wi,i G N, are independent copies of W. The idea of the 
proof is to look at the random path Wi(t) — o~ 2 (t)/2, not from its starting 
point corresponding to t = 0, but rather from its top point. Let us be more 
precise. 

Condition (18) implies that we may define a triple (T, M, F) G R d x R x 
C(R d ) by M = sup tmd (W(t) -a 2 (t)/2), T = ini{t £R d :W(t) - a 2 (t)/2 = 
M} (the "inf" is understood in, e.g., the lexicographic sense) and F(t) = 
W(t + T)-a 2 {t + T) /2 - M. So, (T, M) are the coordinates of the top of the 
path W(t) — a 2 (t)/2, whereas F(t) is the path itself, as viewed from its top. 
Let Mi, Ti and Fi be defined analogously, with W replaced by Wj. Define a 
measurable transformation 

7r : R x C(R d ) -tl^xlx C(R d ) 

by mapping a pair (17, W) G R x C(R d ) to the triple (T, U+M,F) G R d x R x 
C(R d ). Note that E^i fy/i.Wi) is a Poisson point process on R x C(R d ) with 
intensity e~ y dy x dF\y, where Pjy is the law of W on C(R d ). Therefore, by 
the mapping theorem for Poisson point processes (see, e.g., [19]), we obtain 
that 2~2i^=i 3(Ti,Ui+Mi,Fi) is a Poisson point process on R d x R x C(M. d ) whose 
intensity measure \P is given by 

(19) 9(A) = f e- y dyxdP w = [ e- y F[(T,M + y,F)eA]dy, 

Jw- 1 (A) JR 

where A denotes a Borel subset of R d x R x C(R d ). 

We claim that the measure 9 has natural invariance properties. First, it 
follows from (19) that for every z G R, we have 

tf(A + (0,z,0)) = / e- y F[(T,M + y,F)eA + (0,z,0)}dy 

JR 

= f e - y F[(T,M + (y-z),F)eA]dy 
Jr 

= f e~ iy+z) F[{T,M + y,F) €A]dy 
Jr 
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Second, Theorem 10 and Proposition 5 imply that ^(A + (t, 0, 0)) = ^(A) for 
every t £ R d . To see this, note that the collection {(Ti,Ui + M i: Fi),i £ N} can 
be obtained from the collection {Ui + Wj(-) — <r 2 (-), i £ N}, viewed as a trans- 
lation invariant Poisson point process on C(R rf ), by a measurable transfor- 
mation, which commutes with spatial translations. Furthermore, note that 

f([0,lfx[0,l]xC(l d ))< / e~ y F\ sup W(t) > -y] dy 

is finite by the same argument (based on [20]) as in the proof of Proposition 
13. 

We now show that the above invariance properties imply a product- 
type representation for Take a measurable set A C C(R rf ) and con- 
sider a measure ^a on R rf x R, defined as follows: for B C R rf x R, we set 
*&a{B) = J BxA e y d^(t,y, F). By the above, the measure ^ a is translation 
invariant and Vl/^QO, l] d x [0,1]) < oo. It follows that $ a is a multiple of 
the Lebesgue measure and hence we may write dfy a = Q(A)dtdy for some 
finite constant Q(A). Further, Q(A) defines a finite measure on C(R rf ). 
Introducing the normalized measure Q' = Q/c, where c = Q(C(R d )), we may 
write d^f in the form ce~ y dtdy x dQ' . 

We are ready to finish the proof. The Brown-Resnick process of Theorem 
2 may be written as 

oo oo 

V(t) = \J(U t + Wi(t) - a\t)/2) = \/ (F*(t - t*) + y*), 

i=l i=l 

where F*(-) = Fj(-) + logc, t* = Ti and y* = Ui + Mi~ logc. We claim that 
this gives the required mixed moving maxima representation of r\. First, re- 
call that Y^£=\ $(Ti,Ui+Mi,Fi) is a Poisson point process on M. d x R x C(R d ) 
with intensity = ce~ y dtdy x dQ'. It follows that J2t=i $(t*,y*,F*) is a Pois- 
son point process on the same space with intensity e~ y dtdy x dQ* , where 
Q* is the law of F + logc for F ~ Q'. Thus, Y^t=i $(t*,y*) is a Poisson point 
process on R rf x R with intensity e~ y dtdy, whereas F* may be viewed as a 
random mark sampled independently of (t* , y* ) according to the probability 
measure Q*, as required. □ 

Remark 15. In the case d = 1, it follows from Corollary 2.4 of [23] that 
condition (18) is satisfied whenever liminff_ >00 7(t)/logi > 8. 

6. Maxima of independent Gaussian processes. It was shown by Brown 
and Resnick [4] that a suitably normalized and spatially rescaled maximum 
of n independent Brownian motions or Ornstein-Uhlenbeck processes con- 
verges, as n — > oo, to the process 77 of Theorem 1. Some related results were 
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obtained in [15, 17, 18, 24]. We are going to extend the result of [4] to Gaus- 
sian processes whose covariance function satisfies a natural regular variation 
condition. 



Assumption 16. Let {X(t),t G D} be a zero-mean, unit-variance Gaus- 
sian process defined on a neighborhood D C M. d of and having covariance 
function C{t\,t2) = M[X(ti)X(t,2)]- We assume that the asymptotic relation 



(20) 



1 - C{et 1 ,et 2 ) , , 



holds uniformly in £i, £2 G as long as £i , £2 stay bounded. Here, L is a 
function varying slowly at 0, a G (0,2], and 7:M d — > [0,oo) is a continuous 
function satisfying 7(A£) = A Q 7(£) for every A > 0, t G M d . 



Define normalizing sequences 

(21) 6 n = (21ogn) 1 /2_(21ogn)- 1 / 2 ((l/2)loglogn + log(2 v ^)), 

(22) s n =inf{s>0:L(s)s a = 6- 2 } 

and recall (see, e.g., Theorem 1.5.3 in [21]) that, for i.i.d. standard Gaussian 
{Zi,i G N}, we have 



(23) 



lim P 

n— >oo 



\/ -b n ) <y 



i=i 



exp(— e" 



We write r/ n =^ 77 as n — > 00 if, for every compact set K C Mr, the sequence 
of stochastic processes rj n converges to r\ weakly on C(K), the space of 
continuous functions on K. 



Theorem 17. Let Aj,i GN, be independent sample- continuous copies 
of X , a process satisfying Assumption 16. Define 

n 

Vn(t) = \J b n (Xi{s n t) - b n ). 
i=l 

Then, rj n => n as n — > 00, where n is the Brown-Resnick process associated 
to the variogram 27. In particular, 7 must be a variogram. 



Remark 18. The results of [4] can be recovered by applying the above 
theorem to the Ornstein-Uhlenbeck process and to the process X(t) = 
B(t + t)/(t + i) 1 / 2 , where t > and {B(t),t G M.} is a standard Brow- 
nian motion. 
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Proof of Theorem 17. Note that s„,^0asjn oo. Define a process 

Y n {t) = b n {X(s n t)-b n ), tes- l D. 

Further, for losM, let Y™ be the process Y n — w conditioned on {1^(0) = w}. 
Let Yi tn and Y™ n be defined analogously, with X replaced by A^. 

The expectation and covariance of the Gaussian process Y™ are given by 

(24) ^(t) = -(b 2 n + w)(l-C(s n t,0)), 

(25) r n (t 1 ,t 2 ) =b 2 n {C{s n t 1 ,s n t 2 ) - C(s n t 1 ,0)C(s n t 2 ,0)). 

Note that the conditional covariance r n (ti,t2) does not depend on w. Let 
Mi,*2 GR^wGKbe fixed. Using (20) and (22), we obtain 

(26) Hm tf(t) =- 7 (t), 

n — »oo 

(27) lim r fl (t 1 ,« 2 )=7(ti)+7(t 2 )-7(*i-*2). 

n— »oo 

A further consequence of (24) is that as long as t stays bounded, there is a 
constant c such that, for sufficiently large n, we have 

(28) \n™(t)\<c+\w\/2 VroGR. 

It follows from (26), (27) that as n— > oo, the process Y™ converges in 
the sense of finite-dimensional distributions to {W(£) — j(t),t £ M. d }, where 
{W(t),t £ M d } is a Gaussian process with stationary increments, variogram 
27 and TU(0) = 0. On the other hand, it is known (see, e.g., Corollary 4.19 
in [26]), that the point process Ya=i „(o) converges, as n — > 00, to the 
Poisson point process on M. with intensity e~ y dy. From these two facts, at 
least on the formal level, we obtain the statement of the theorem. However, 
making this rigorous requires some work. 

First, we show that ij n converges to r] in the sense of finite-dimensional 
distributions. Let t\,...,tk 6 R rf and yi,.-.,yk SRbe fixed. By conditioning 
on Y n (0) = w and noting that the density of Y n (0) is given by 

/y„(o)H = l/(^n)e-^ +b ") 2 / 2b ", 

we obtain 

Ppj:Y n (t j )>y j ] 

e -^+bl) 2 /(2bl) F ^ j . Yn ( tj ) > y j \Y n (0)=w} dw 



1 



27r6 n 
1 



2irb n e b2 J 2 



^—w—w 1 



/(2&£)p[ 3j - . Y™{tj) > Vj - w] dw. 
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Noting that (21) implies that y/2irb n e b ™/ 2 ~ n as n — > oo and taking A > 0, 
we may write the above as 

P[3j:Y n (t j )>y J ]r,lff A + r + f = h h {n) + h{n) + I 3 (n)). 
n \J-A J A J-oo J n 

Since the convergence of the distribution of {Y™{ yj )} k j=1 to that of {W(yj) - 
l(yj)}j=i 1S uniform provided that w G [— A, A], we obtain 

rA 

(29) lim Ji(n) = / e~ w F[3j : W(tj) - */(tj) > yj - w] dw. 

n^oo J— A 

For l2(n), we have the trivial estimate 

;>oo 

(30) h{n) < / e~ w dw = e- A . 



We estimate Iz{n). Using (28), we obtain, if w < —A and A,n are large, 
nYn(tj) > Vj -w}< F[Y^( tj ) - > yj -c- M/2 - w] 

tf(t j )>\w\/4]. 

Recall the well-known estimate ^f(t) < e~* / 2 , t > 0, where ^(t) is the tail 
of the standard Gaussian distribution. By (27), Var[Y^" (?;,■)] < k 2 for some 
K > and all j = 1, . . . , k, n G N. Hence, 

p[y n w (ty) > -w] < e -(™/ 4 ) 2 /(2« 2 ). 

It follows that 

ia(n)<5Z / A e- w F[Y^(t j )>y j -w]dw<k f * e~ w e~ w * I V 2 "^ dw . 

j_^J—00 J— oo 

Hence, 

(31) lim limsup 13(71) = 0. 

A— >oo n^oo 

Bringing (29), (30) and (31) together and letting A — > 00, we obtain 
F[3j : Y n (tj) > Vj ] ~ - / e - w P[3j : Wfe) - 7^) > yj - w] 



n 

= —Eexp max (W(tj) — j(tj) — yj) 
n j=i,...,k 

as n — > 00. Therefore, 

lim P[Vj : rinitj) < y 3 ] = lim (1 - F[3j : Y n (t 3 ) > 



exp< —Eexp max (W(tj) — j(tj) — yj) \. 
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By (5), the right-hand side coincides with :rj(tj) < yA, which proves 
that r) n converges to r/ in the sense of finite-dimensional distributions. 

It remains to show that the sequence rj n is tight in C(K), where K C M. d 
is a fixed compact set. First, note that the sequence r] n (0) is tight in R [in 
fact, the distribution of ry n (0) converges weakly to the Gumbel distribution]. 
For a function / G C{K) and 5 > 0, define 

cj 5 (f) = sup |/( tl )_/( t2 )|. 

tl,t2GK,||ti-t 2 ||<5 

By the standard tightness criterion (see, e.g., Theorem 7.3 in [3]), we need 
to show that for every e > 0, a > 0, there exists some 5 > such that 

(32) P[<4jfan) >a]<e for all n > N. 

Throughout, N denotes a large integer whose value may change from line 
to line. We concentrate on proving (32). The proof of the next lemma will 
be given later. 

Lemma 19. The following assertions hold: 

1. for every c > 0, the family of processes Y™ , w G [— c,c], n£N, is tight in 
C(K); ' 

2. the family of processes Y™ — [i™, w£M, n G N, is ii<?/ti in C(K). 



For ci > 0, define a sequence of random events 

E n = { inf 77 n (t) < -ci}. 

We show that we can find c\ > such that P[£? n ] < e for all n> N. First, 
choose Co so large that 2e _c ° < e. Using part 1 of Lemma 19, choose c\ so 
large that 

mTY™(t) < Co - ci <l/2 for all w G [-c , c ], n G N. 
Define random events 

A,n = |V;,n(0) G [-C0,0o], ilrf Fi, n (t) - li,n(0) > C - Ci}. 

We have, by conditioning on 5^ jn (0) = w, 



P[Ai, n ] = (V^bne^y 1 f C0 e~ u - w 2 /( 26 «)p[inf Y™(t) > c - d 



4n ./-«, 
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which implies that P[A$ n ] > cq/u if cq is sufficiently large and n > N 
N(c ). Noting that P[£ n ] < P[n"=i A? n ] gives 



P[£„] < (1 - co/n) n < 2e~ c ° < e, 
For C2 > 0, define the random events 



n> N. 



F n = \j{Yi,n(0) >C 2 }, 
i=l 

G n = {3teK:7 ]n (t)^ 



sup *i,n(*)f- 
ie{l,...,n}: |y i>n (0)|<C2 J 



Trivially, P[-F n ] = P[r/ n (0) > c 2 ] < e for every n, if c 2 is large. We show that 
there exists some c 2 > such that P[G n ] < 3e for n > N. Introduce random 
events 

Bi,n= !^,n(0) < -C 2 ,SUpYj n (t) > -Cl\. 

L teK ' 

Then, again conditioning on Yi tTl (0) = w and recalling (28), we obtain 

P[B it „] = (V^bne^ 2 )- 1 [ ^ e-^-^/^P [sup Y^i) > -c x - w] dw. 

J-00 HeK J 

By part 2 of Lemma 19, there exists some C3 > such that 

sup(Y™(t) -AC(i)) >c 3 l <l/2, wel,nGN. 
teAT J 

Recall that, by (28) and (25), we have 

sup//™<c 4 , supVarY™(i) < k 2 , w<0,n>N, 

teK 2 jgA- 

for some C4,k. Applying Borell's inequality (see Theorem D.l in [25]), to- 
gether with the above estimates, we obtain, for w < 0, 

sup Y™(i) >-a-w < 2*(-(-ci - w/2 - c 3 - c 4 )/k), 
teK J 

where \I/ is the tail of the standard Gaussian distribution. If w < — 4(ci + 
C3 + C4), this, together with the bound ^f(t) < e~ l / 2 , t > 0, implies that 



supY^(i) > -ci - u; 
teA 



< 



2e - W V(32 K 2 ) n>Ar> 



It follows that, for n> N and c 2 > 4(ci + C3 + C4) 



[Si,„] < 



7? 



-C2 



e e 
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So, we can choose c 2 sufficiently large that nP[£>i in ] < e for n> N . Therefore, 
P[G„] < F[E n ] + F[F n ] + P[G n \ (E n U F n )\ <2e + n¥[B hn ] < 3e. 
We are now ready to prove (32). Let 

Ci, n = {*i,n(0) G [-C2,c 2 ],CJ5(yj in ) > a} 
and define H n = UILi Ci,n- Then, 

P[Q, n ] = (v^fone 6 '/ 2 )- 1 P e-— 2 ^ 26 "^^^^) > a] dtc. 

J— c 2 

By part 1 of Lemma 19 and the tightness criterion (see Theorem 7.3 in 
[3]), we can make F[ug(Y™) > a] arbitrary small (uniformly in w G [—02,02] 
and for n > N) by choosing 5 small. So, choose 5 > sufficiently small that 
F[C i>n ] < f . Then, 

F[iv s {Vn) >a}< F[G n ] + F[H n ] < 3e + nP[G\ n ] < 4e, 

which yields (32) with 4e instead of e. This proves the tightness of the 
sequence rj n and completes the proof of Theorem 17. □ 

PROOF of Lemma 19. It follows from (25) that, independently of w, 

VarCTO-Y-fe)) 

= b 2 n (2 - 2C(s n t 1 ,s n t 2 ) - (C(s n ti,0) - C(s n t 2 ,0)) 2 ) 
<2b 2 n (l-C(s n t 1 ,s n t 2 )). 

Assumption 16 implies that, uniformly in t\,t 2 G K, 

Vax(l^(ti) - y n w fe)) < 26 2 • 2(L( 8n )a« 7 (i 1 - t 2 )), n>N. 

By (22), we have b 2 l L(s n )s° <2,n>N, and so, for some C5 > 0, 

(33) Var(y r r(t 1 )-y n w (t2))<87(ti-t2)<c 5 ||ti-t2|| Q , n > N. 

Now, the second claim of the lemma follows from (33) by applying Corol- 
lary 11.7 of [22] to the family of processes Y™ - /i™, w £R, n G N [take 
-0(x) = x 2 , d 2 {t\,t 2 ) = cspi — £2|| Q there]. To prove the first claim, we need 
to additionally show that [j,™, w G [— c, c], n G N, is a tight family of func- 
tions in C{K). This last statement follows from (24), which shows that the 
convergence fJ%(t) — > — "f(t) in (26) is uniform in t G K, w G [— c,c]. □ 
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7. Domains of attraction. We are now going to prove a partial con- 
verse of Theorem 17. More precisely, we characterize all nontrivial limits 
of normalized and spatially rescaled pointwise maxima of stationary Gaus- 
sian processes. Let us call a random process {rj(t),t degenerate if, for 
all t\,t<i G we have r}{t\) = rjfo) a.s. 

Theorem 20. Let {X(t),t G be a stationary zero-mean, unit-variance 
Gaussian process with continuous covariance C(t) = K[X(Q)X(t)] and let Xi, 
iEN, be independent copies of X . Suppose that, for some sequences a' n > 0, 
b' n G M. and s' n > 0, the process {n' n (t),t G M d } defined by 

n 

Vn(t) = ya n (X t (s' n t)-b' n ) 

i=l 

converges, asn^oo, to some nondegenerate, continuous-in-probability pro- 



cess {rj'(t),t G R }, in the sense of finite- dimensional distributions. The fol- 
lowing assertions then hold: 

1. there is an a G (0,2], a finite measure \i on the unit sphere S^" 1 in M. d 
and a function L that varies slowly at such that 

(34) l-C(t)~L(||i||) 7 (i) ast^O, 

where 



(35) 7 (t)=/ |(t,x)| Q d^); 

2. £/te normalizing sequences a' n ,b' n and s' n satisfy 

(36) lim a' /&„ = A > 0, lim b n (b' - b n ) = B € R, 

n— »oo n— >oo 

(37) lim &^(4)< Q = s > 0, 



where b n is defined by (21); 
3. the limiting process r)' coincides with A(t] — B), where rj is the Brown- 
Resnick process associated to the variogram 2s"f. 

We need a lemma, the essential part of which was proven in [18]. 

(n) (n) 

Lemma 21. For n G N, let Z\ ,...,Z n be i.i.d. bivariate Gaussian 
vectors having standard Gaussian margins and correlation p n . The maxima 

n 

M n =\/ b Tl (Z\ n) -b n ) 
i=l 

converge in distribution to some bivariate random vector if and only if 
(38) lim b 2 n (l -p n ) = c 
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for some c£ [0,oo]. The limiting bivariate distribution depends on c contin- 
uously; its margins are independent if and only if c = oo and are equal a.s. 
if and only if c = 0. 

Proof. Suppose, first, that (38) holds. Then, by a result of [18], the 
sequence M n converges in distribution. The explicit formula, given in [18], 
shows that the limiting distributions corresponding to different values of 
c are different. Suppose, now, that (38) does not hold. We then have < 
liminf 6^(1 — p n ) < limsup6^(l — p n ) < oo. Again using [18], we obtain that 
the sequence M n has at least two different accumulation points and thus does 
not converge. The last claim of the lemma follows again from the explicit 
formula in [18]. □ 

Proof of Theorem 20. By stationarity of X, the distribution of rf{t) 
does not depend on t G M. d . Thus, if for some constant Co, f/(0) = Co a.s., then 
for every t € W 1 , rf(t) = Co a.s., which is a contradiction since rf is assumed 
to be nondegenerate. So, in the sequel, we assume that r/(0) is not a.s. 
constant. In this case, the convergence-to-types theorem (see Proposition 
0.2 in [26]), together with (23), yields constants 4>0,fi£t such that (36) 
holds. It follows that the process 

n 

Vn(t) = V b n (Xi(s' n t) - b n ) 
i=i 

converges, as n — ► oo, to the nondegenerate limit rj = A" 1 ?]' + B. From now 
on, we consider the processes r] n and rj instead of rj' n and rj'. 

For any fixed t £ M. d , the previous lemma, applied to the triangular ar- 
ray of bivariate vectors z\ n ^ = pQ(0), Xi(s' n t)), i = 1, . . . ,n, n € N, yields a 
constant c(t) G [0, oo] such that 

(39) \mijl{l-C{s l n t)) = c{t). 

Since the limiting process rj is assumed to be continuous in probability, the 
distribution of the bivariate vector (t](0), Tj(t)) must converge weakly to the 
distribution of (77(0), 77(0)) as t — > 0. Using the last statement of Lemma 21, 
we obtain that limj_>o c(t) = c(0) = 0, that is, c is continuous at the origin. 
Note, also, that c(io) 7^ for some to 7^ since otherwise the process r/ would 
be degenerate. 

By Bochner's theorem, there exists an IR^-valued random variable £ such 
that the characteristic function of £ is C(t). Moreover, since the function C 
is real-valued, the distribution of £ must be symmetric with respect to the 
origin. Let £j, % G N, be i.i.d. copies of £. Then, the characteristic function 

ip n Of 

K] 
i=l 
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is given by <p n (i) = C(s' n t)^ so that (39) yields 

lim y n (t) = lim (1 - c(t)/b 2 n + o{l/b 2 n )ftt = e -(*). 

n — >oo n — >oo 

Now, Levy's convergence theorem tells us that the random vector S n 
converges weakly to a random vector 5 whose distribution is necessarily 
nondegenerate (i.e., F[S = 0] / 1; to see this, recall that c(to) ^ and hence 
g-c(to) ^ i f or som e to 7^ 0), stable with some parameter a € (0,2] and sym- 
metric with respect to the origin. It follows from the characterization of 
domains of attraction of multidimensional symmetric stable distributions in 
terms of characteristic functions (see Corollaries 1 and 2 in [1]) that the 
covariance function C must have the form (34), (35). Inserting this in (39) 
for some t with = 1, we obtain 

fim b 2 n L(s' n )s'« 1 (t) = c(t), 

which yields (37). Furthermore, (34) and (35) imply that the process X sat- 
isfies Assumption 16. Therefore, by Theorem 17, the limiting process rj must 
be the Brown-Resnick process associated to the variogram 257. Recalling 
that i] = A~^r( + B, we obtain the last statement of the theorem. □ 

8. Extensions and remarks. In view of Theorems 17 and 20, the ques- 
tion arises as to whether max-stable processes corresponding to variograms 
7 that are not of the form (35) also admit representations as limits of point- 
wise maxima of stationary Gaussian processes in some broader sense, as in 
Theorem 20. The answer is affirmative, as the following theorem shows. 

Theorem 22. Let 7 be a variogram on M. d , that is, 7(0) = and 7 is 
negative definite. For each n G N, let X\ n , . . . ,X nn be i.i.d. copies of a sta- 
tionary zero-mean Gaussian process {X n (t),t £ with covariance function 
exp(-7(i)/f£). Define 

n 

Vn(t) = \/ b n (X m (t)-b n ), t£R d . 

i=l 

Then, rj n converges, in the sense of finite- dimensional distributions, to the 
Brown-Resnick process associated to the variogram 27. 

Proof. Note that exp(— 'y(t)/bf l ) is indeed a covariance function of some 
stationary Gaussian process, by Schoenberg's theorem (see Theorem 7.8 in 
[2]). As in the proof of Theorem 17, it can be shown that the conditional 
distribution of b n (Xi n (t) — Aj„(0)), given that b n (Xi n (0) — b n ) = w, converges 
to the distribution of W(t) — j(t), where If is a Gaussian process with 
stationary increments, variogram 27 and Vf(0) = 0. The rest of the proof is 
the same as that of Theorem 17. □ 
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Remark 23. The above theorem gives another proof of stationarity in 
Theorem 2 in the case W(0) = 0. 

Remark 24. The bivariate distributions of the Brown-Resnick process 
T) associated to the variogram 7 are given by the formula 



where <3? is the standard normal distribution function. 

Proof. The remark is a consequence of Theorem 22 and a result of 
[18]. Moreover, it follows from Theorem 22 that the finite-dimensional dis- 
tributions of the process r\ belong to the family of multivariate max-stable 
distributions introduced in [18]. □ 

Remark 25. A natural dependence measure between 7/(0) and rf(t) is 
given by p(t) = 2 — s(t) E [0, 1], where <j(t) is determined from the condition 



for some (and hence all) zfl; see, for example, [7, 28]. It follows from 
Remark 24 that 



Thus, a variogram 7 is completely determined by the dependence function 
p(t) of the corresponding process rj. It follows that 7/(0) and rj(t) become 
asymptotically independent as ||t|| — > 00 [which corresponds to p(t) — > 0] if 
and only if y(t) — > 00 as ||t|| — > 00. Furthermore, if d = 1, then, by Theorem 
3.4 in [30], the process 77 is mixing if and only if 7(4) — ► 00 as t — > 00. 

Remark 26. Theorem 17 may be generalized to processes whose co- 
variance has different Holder exponents in different directions. For example, 
assume that {X(t),t G M. d } is a stationary zero-mean Gaussian process with 
covariance function C satisfying 



nr ] (t 1 )<y 1 ,r ] (t 2 )<y 2 ) 





P[t?(0) < z, r](t) <z]= Pfa(0) < z] q{t) 



p(i) = 2(l-$(^)/2)). 



C(t) = C(ti, . . . ,t d ) = 1 - ^ c t \ur + o(\\t\\ ad ) 



as i 







i=i 
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for some < a\ < ■ ■ ■ < a d < 2, c\, . . . , c<j > 0. If Xi,i € N, are independent 
copies of X, then 

n 

Vn(t) = V bniXiib-^h, b- 2 / a H d ) - b n ) 
i=l 

converges to the Brown-Resnick process associated to the variogram 2j, 
where 7(ii, . . . ,t d ) = £f =1 Ci\U\ ai . 
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